On capillary-gravity waves generated by a slow moving object 
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We investigate theoretically and experimentally the capillary-gravity waves created by a small 
object moving steadily at the water-air interface along a circular trajectory. It is well established 
that, for straight uniform motion, no steady waves appear at velocities below the minimum phase 
velocity c m i n = 23 cm • s _1 . We show theoretically that no such velocity threshold exists for a steady 
circular motion, for which, even for small velocities, a finite wave drag is experienced by the object. 
This wave drag originates from the emission of a spiral-like wave pattern. Our results are in good 
agreement with direct experimental observations of the wave pattern created by a circularly moving 
needle in contact with water. Our study leads to new insights into the problem of animal locomotion 
at the water-air interface. 

PACS numbers: 47.35.-i , 68.03.-g 



Capillary-gravity waves propagating at the free surface 
of a liquid are driven by a balance between the liquid in- 
ertia and its tendency, under the action of gravity and 
surface tension forces, to return to a state of stable equi- 
librium For an inviscid liquid of infinite depth, the 
dispersion relation relating the angular frequency lo to 
the wave number k is given by u> 2 = gk + jk 3 / p, where p 
is the liquid density, 7 the liquid-air surface tension, and 
g the acceleration due to gravity • The above equation 

may also be written as a dependence of the wave velocity 

1/2 

c(k) = ui{k)/k on wave number: c(k) = (g/k + jk/ p) ' . 
The dispersive nature of capillary-gravity waves is re- 
sponsible for the complicated wave pattern generated at 
the free surface of a still liquid by a moving disturbance 
such as a partially immersed object {e.g. a boat or an in- 
sect) or an external surface pressure source [1, H, 0, HI, El • 
Since the disturbance expends a power to generate these 
waves, it will experience a drag, R w , called the wave re- 
sistance [3]. In the case of boats and large ships, this 
drag is known to be a major source of resistance and im- 
portant efforts have been devoted to the design of hulls 
minimizing it Q- The case of objects small relative to 
the capillary length k _1 = (ll(pg)) 1 ^ 2 has only recently 
been considered 

In the case of a disturbance moving at constant veloc- 
ity V, the wave resistance R w cancels out for V < c min 
where V stands for the magnitude of the velocity, and 
{ig-y/p) 1 / 4 is the minimum of the wave velocity 
c(k) given above for capillarity gravity waves [!, 0, Hf. 
For water with 7 = 73 mN • m -1 and p — 10 3 kg ■ m~ 3 , 
one has c mul — 0.23 m • s -1 (room temperature). This 
striking behavior of Rw around c m ; n is similar to the 
well-known Cerenkov radiation emitted by a charged par- 
ticle and has been recently studied experimentally 
13, 14]. In this letter, we demonstrate that just like ac- 
celerated charged particles radiate electromagnetic waves 
even while moving slower than the speed of light [l5l ]. 
an accelerated disturbance experiences a non-zero wave 



resistance R w even when propagating below c m i n . We 
consider the special case of a uniform circular trajec- 
tory, a situation of particular importance for the study 
of whirligig beetles (Gyrinidae, [161 ]) whose characteris- 
tic circular motion might facilitate the emission of surface 
waves that they are thought to be used for echolocation 
[13, EH ■ This work is therefore restricted to the effect of 
a wake stationary in the rotating frame, and do not con- 
sider time dependent contributions, like vortex shedding 

HI S3. 

We consider the case of an incompressible infinitely 
deep liquid whose free surface is unlimited. In the ab- 
sence of external perturbation, the free surface is flat 
and each of its points can be described by a radius vec- 
tor r = (x, y) in the horizontal plane. The motion of a 
small object along the free surface disturbs the equilib- 
rium position of the fluid, and each point of the free sur- 
face acquires a finite vertical displacement £( r )- Rather 
than solving the complex hydrodynamic problem of find- 
ing the flow around a moving object, we consider the dis- 
placement of an external pressure source P ext (r, t) H,[l]. 
The equations of motion can then be linearized in the 
limit of small wave amplitudes [l9| . 

In the frame of this linear-response theory, it is conve- 
nient to introduce the Fourier transforms of the pressure 
source P ext (k,t) and of the vertical displacement ((k,t) 
[2fj| ] . It can be shown that, in the limit of small kinematic 
viscosity v, the relation between C(fc,0 an d P ex t(k,t) is 
given by 



d 2 C , A , 2 d( , a , M ; kP ext (k,t) 

W + Avk at + u (fc)c = p — 



(i) 



In this letter we assume that the pressure source has 
radial symmetry and that the trajectory ro(i) of the ob- 
ject is circular, namely : r (t) = 1Z (cos(Qi), sin(f2t)). 
Here 7Z is the circle radius, and ft is the angular fre- 
quency. The linear velocity of the object is then given by 



V = TZCl. With these assumptions, the external pressure 
field isP ext (r,t) = P ext (\r - r (t)\, yielding in Fourier 
space P ext (k,t) = P ext (fc)e _lfc r ° (t) . Since the right hand 
side of Eq. ^ is periodic with frequency SI, it is possible 
to find its steady state solution by expanding the right 
hand side into Fourier series. The problem then becomes 
equivalent to the response of a damped oscillator to a 
sum of periodic forces with frequencies nil, where n is 
an integer. The vertical deformation at any time t can 
then be reconstructed by evaluating the inverse Fourier 
transform. For the particular case of uniform circular 
motion, the time dependence is rather simple. Indeed, 
in steady state, the deformation profile rotates with the 
same frequency f2 as the disturbance. Therefore, in the 
rotating frame, C depends on the position r only. The 
analytical expression of C( r ) m cylindrical coordinates 
(x, y) = r(cos (/>, sin cj>) is given by 



k 2 dk P ext (k)J n (kr)J n (kll) 
2np n 2 Vl 2 - u 2 (k) + Ainvk 2 Vl 



(2) 



where J n is u-th order Bessel function of the first kind. 
The summation index n is directly related to the n-th 
Fourier harmonic of the periodic function e -* fc r o(*) anc ^ 
since the problem is linear, the contributions of all the 
harmonics add together. 

The knowledge of the exact structure of the wave pat- 
tern is precious, but a quantitative measurement of the 
wave resistance is needed in order to understand, for ex- 
ample, the forces developed by small animals moving at 
the surface of water. In the case of the circular motion 
under study, the wave resistance R w can be calculated 

from its average power P w — — J d 2 r (^P ex t {r, t) 9 ^'^ ^ 

by R w — P w /V. Using the Fourier expansion of (, 
one then obtains in the limit ^At/c m i n — > (for water, 
w/cmi n ~ 10" 3 ): 



R W (V,TZ) 
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n>0 
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(3) 



dk Jk„ 



where k n is the unique solution of the equation oj{k n ) = 
nfi (the notation R W (V,1Z) stresses the dependence of 
R w on the velocity magnitude and on the trajectory ra- 
dius). Equation ([3]) shows that the wave resistance R w 
takes the form of a sum R w = X^n>o where the A n 
are positive numbers that measure the contribution of 
each Fourier mode of the external pressure source (with 
frequency n il) to the wave resistance. 

A numerical calculation of the wave resistance is pre- 
sented in Fig. [T] for a pressure source P ext {k) — 
Poexp(-kb), where po is the total force exerted on the 
surface and b is the typical object size (22j. As observed, 
Eq. [3] differs significantly from the original prediction on 
the wave drag in the case of a straight uniform motion 
with velocity V [1, [ll| given by 




V/c n 



FIG. 1: (Color online) Plot of the wave resistance R w in 
units of Poii/"f, as a function the reduced velocity V/c m m = 
7?.f2/c m i n for different ratios between the trajectory radius 1Z, 
and the object size b, as predicted by Eq. (|3}- The red curve 
(presenting many oscillations) corresponds to IZ/b — 100, 
while the black one (with fewer oscillations) corresponds to 
IZ/b = 10. The green curve displaying a typical discontinuity 
at V = c m jn is the wave drag for a straight uniform motion 
with velocity V [!|. The object size, b, was set to b = 0.1 ft" 1 . 



Rw,i(Y) = 



kdkP 2 xt {k)9{V 



z{k)) 



o 2tt P V 2 y/l-(c(k)/Vy 



(4) 



where 9(.) is the Heavyside function and c(k) = u>(k)/k 
is the phase velocity. Most notably, the wave drag 
for a circular motion is non-zero for all velocities, even 
for V < c m j„ where wave-resistance vanishes exactly in 
the case of a linear motion and this effect is far from 
negligible: for IZ/b = 10 and at velocities as slow as 
V/c m i n ~ 0.6, the wave drag is still one fifth of that 
applied to an object moving linearly at V/c m i n = 1. 
The radiation of waves by an accelerated particle should 
not be surprising and actually is a very general phe- 
nomenon that can be observed for instance in elec- 
tromagnetism (bremsstrahlun q) o r in general relativity 
(Zeldovich-Starobinsky effect [23|]). Mathematically, the 
fact that, for a circular motion, the wave resistance is 
finite even below c m - m can be understood as follows. In 
the case of uniform motion, all the wavenumbers such 
as c(k) < V contribute to the wave drag, whereas for 
circular motion this is the case for only a discrete set 
of wavenumbers fc„. While the condition c(k) < V can 
be satisfied only when V > c m ; n , the equations for the 
wavenumber k n , u>(k n ) = riV/lZ, have positive solutions 
for any velocity V . These wavenumbers k n create finite 
contributions A n > to the wave drag. Therefore for a 
circular trajectory a finite wave drag exists at any veloc- 
ity V > 0; for the same reasons R w is also continuous at 
V = c m j n . Moreover, the wave resistance develops a small 
oscillating component as a function of the velocity V. It 
originates from the oscillatory behavior of Bessel func- 
tions and will be analyzed more thoroughly in a future 
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FIG. 2: (color online) Wave radiation for V « 21cm /s 
ss 0.9c m in with a radius 1Z « 2.7cm ~ 9k _1 Left: Color di- 
agram of the surface deformation £(r) computed numerically 
from Eq. @. This image represents a square region of size 
400re~ around the center of rotation, red color corresponds 
to maximal £(r) values, while green corresponds to minimal 
values of £(r). The cross indicates the center of the trajec- 
tory and the moving object is located in the region of highest 
deformation. Right: Photography of the wave crests gener- 
ated on a water surface by a needle rotating at a velocity . On 
both pictures, the black curve represents the the Archimedean 
spiral of radius given by Eq. (0 



publication. Finally, we note that despite these strik- 
ing differences Eqn. |3]) and Q should coincide in the 
limit of a large trajectory radius TZ. We confirmed this 
behavior by checking both analytically [24| and numer- 
ically that in the limit TZ — > oo, R W (V,1Z) — > R w ,i(V). 
However even if the circular wave drag R W (V, TZ) is close 
to R Wt i{V) starting from TZ/b ~ 10, important differ- 
ences remain even up to very large values of TZ/b such as 
TZ/b ~ 100. 

Figure 2 represents the wave crest pattern (computed 
numerically form Eq.Q) at the origin of this finite wave 
drag. It exhibits characteristic concentric Archimedean 
spirals (also known as arithmetic spirals) of the form 
r = a0 + ro. This can be understood from our theoretical 
results as follows. In a first estimation, one can assume 
that the integrals in equation Eq. ^ are dominated by 
the contribution of the poles at k = k n . Thus £(r) can 
be written as £(r) ~ A= J2 n B n e l ^~ knr \ where we have 
used the asymptotic development of J n {k n r) at large dis- 
tances r and B n are complex coefficients that do not de- 
pend on the position r = r(cos</>, sin<^>). By separating 
the contribution of the different modes in the relation 
F(t) = -J d 2 rP ext (r,t)V((r,t), one finds that B„ is 
proportional to A n (where, as defined earlier, the positive 
coefficients A n measure the contribution of each Fourier 
mode to the wave drag: R w = J2 n >o A n )- One can show 
that in the regime of small object sizes del, the pro- 
portionality constant between B n and A n depends only 
weakly on the Fourier mode number n; thus, one has 
C(r) oc A= J2 n A„e 4 '™^ _fc " r '. We have checked numeri- 
cally that in the regime V < c m i n , the distribution of the 
coefficients A n is usually peaked around n ~ kTZ. For ex- 



ample, for kTZ = 10 and nb = 0.1, A n is peaked around 
n = 10 for velocities V in the interval (c min /2, c min ). 
The wave-crests are given by the lines of constant phase 
n<p — k„ r — const of the dominant mode n = kTZ, leading 
to the following expression for a: 



kTZ 



k{uo = kV) 



(5) 



where k{uj) is the inverse function of u>(k). An interest- 
ing special case of the formula Eq. ([5]) corresponds to 



V 



for which one obtains a m TZ. The spiral 



predicted by Eq. |(5J) is in very good agreement with the 
exact numerical results (Eq. @), as can be seen in Fig. 2. 

We have also compared our theoretical approach with 
experimental results obtained using a one millimeter wide 
stainless steel needle immersed in a 38 cm wide water 
bucket. The needle was rotated on circular trajectories 
of various radii and angular velocities. Since direct mea- 
surement of wave drag, and in particular comparison with 
theory, is non-trivial even for a linear motion [l3l . [l4| , we 
restricted ourselves to the study of the wake itself. A typ- 
ical wave pattern obtained by this method is shown on 
Fig. 3 for TZ ~ 2.7cm and £1 » 2ir x 1.2 Hz (correspond- 
ing to V/c m i n ~ 0.9) and unambiguously demonstrates 
the existence of a wake at velocities smaller than c m i n . 
The observed wave pattern is in remarkable agreement 
with the theoretical prediction r = a<f> + b with a given 
by Eq. ([5]) and r$ a free parameter corresponding to an 
overall rotation of the spiral [25| . For V/c m i n lower than 
0.8, no wake was observed by naked eye. At lower ro- 
tation velocities, we probed the surface deformation by 
measuring the deflexion of a laser beam reflected by the 
air- water interface at a distance r = 11 cm from the ro- 
tation axis. 

Using this scheme, we have established the existence 
of waves down to V/c m i n « 0.6, and verified quantita- 
tively that the wave packet spectrum is peaked around 

< ijj >~ kTZ f2 (see Fig. 3). Experimentally, the fre- 
quency < ui > corresponds to the period of the fast tem- 
poral oscillations of the laser deflection angle (see Fig. 3 
inset). In order to compare our experimental results with 
our model, we note that the deflection of the laser at a 
point r is proportional to the derivatives ~ g|r^ an( ^ 

d ^dr ^ " ^ or siraplicityi we will mainly consider the angu- 
lar derivative, but we have checked numerically that our 
result do not depend on this choice. Using Eq. {2J the 
angular derivative can be decomposed into Fourier se- 
ries: = J2 n C n e in ^- nt ^ lk " r . The coefficients C n 
are proportional to the contribution of the frequency nfl 
to the wave packet spectrum and we can thus calculate 
the mean wave packet frequency using the expression: 

< lj >= ^^2 n>0 n\C n \/^2 n>0 \C n \. As shown in Fig. 3, 
our model is consistent with good accurcy with the ex- 
perimental data without any adjustable parameters. 

Below V/cmin ~ 0.6, the signal to noise ratio of the 
experiment becomes to small to observe the laser deflec- 
tion. Note that this value is in qualitative agreement 
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FIG. 3: (Color online) Inset: typical time dependence of 
the laser deflection angle (arbitrary units) during a rotation 
period T ro t = 2n/Q, the fast oscillation frequency is given 
by < uj >— 2tt/T. Main figure: Dependence of the ratio 
< ui > /fi on kR, for different needle velocities. The dashed 
curves represent experimental results, while the continuous 
curve display the numerical results of our model. Red, green 
and blue curves (diamonds, squares and circles respectively) 
correspond to U/c m i n = 0.69, 0.76 and 0.84. The black curve 
correspond to the analytical estimate < ui > /fi = kTZ. 

with Fig. 1 where the wave resistance (hence the wave 
amplitude) has also significatively decreased with respect 
to its maximum value for V/c m i n < 0.5: we indeed note 
that for 



To summarize, we have shown theoretically that a dis- 
turbance moving along a circular trajectory experienced 
a wave drag even at angular velocities corresponding to 
V < c m j n , where c m i n is the minimum phase velocity 
of capillary-gravity waves. Our prediction is supported 
by experimental observation of a long distance wake for 
V/c m i n as low as 0.6. For V/c min > 0.8, we observed 
by naked eye Archimedean spiral shaped crests, in good 
agreement with theory. These results are directly related 
to the accelerated nature of the circular motion, and 
thus do not contradict the commonly accepted thresh- 
old V = c m i n that is only valid for a rectilinear uniform 
motion, an assumption often overlooked in the literature. 
It would be very interesting to know if whirligig beetles 
can take advantage of such spirals for echolocation pur- 
poses. Although restricted to stationary wakes and thus 
excluded effects such as vortex shedding, the results pre- 
sented in this letter should be important for a better 
understanding of the propulsion of water- walking insects 
[H, [U [H, [29| where accelerated motions frequently oc- 
curs (e.g when hunting a prey or escaping a predator 
[3fj|). Even in the case where the insect motion is rec- 
tilinear and uniform, one has to keep in mind that the 
rapid leg strokes are accelerated and might produce a 
wave drag even below c m ; n . 
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